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Light hadron spectroscopy with O(a) improved dynamical fermions. 

M. Talevi a *, for the UKQCD Collaboration 

a Department of Physics & Astronomy, University of Edinburgh, The King's Buildings, EH9 3JZ (UK) 

qq ' We present results for the hadron spectrum and static quark potential from a simulation with two flavours 

Q\ , of O(a) improved dynamical Wilson fermions at /? = 5.2. We address the issues of sea quark dependence of 
^\ ■ observables and finite-size effects. 
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Quenching is the only non-systematically im- 
provable approximation in lattice QCD. More- 
over, sea-quark effects are to be taken into ac- 
count if we are to make theoretical predictions 
on such phenomena as rf — 7r splitting and string 
breaking. Dynamical simulations are computa- 
tionally very expensive and control of system- 
atic effects is thus of the utmost importance. 
The improvement counterterm csw(<?o) for the 
SW-improved action has been computed non- 
perturbatively, thus yielding full 0{a) improve- 
ment with Nj = 2 dynamical fermions 0. 

The simulation parameters are summarized in 
tab. 0. For each volume and each value of 

^sca 

we have simulated valence quarks with an ap- 
propriate set of k's to exploit the full (K sea , K va i) 
plane for "strange" physics At {3 = 5.2 the 
PS/V mass ratios for the relevant k's are shown 
in tab. ||. The value of csw used differs from 
the final estimate of the Alpha Collaboration 
but since a notable part of the simulation 
had already been carried out with this improve- 
ment coefficient, UKQCD have decided to con- 
tinue the simulation, concentrating on finite-size 
effects and looking for signs of sea-quark effects 
in both the static quark potential and the light 
hadron spectrum. The complete removal of 0(a) 
artefacts is thus not the highest priority in this 
work. Details of the simulation and results can 
be found in [Q], while we refer the reader to || 
for an up-to-date review on spectroscopy. 
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Simulation parameters. 

1. Static quark potential 

To set the lattice spacing, we have used the 
scale ro determined from 

F(r /a)(r /a) 2 = 1.65, r = 0.49 fm, (1) 

where F(r) is the force between two static quarks, 
derived from the potential V(r), which is in turn 
extracted in a standard way from Wilson loops 
H . In Full QCD r is better suited than the string 
tension, since the string is expected to break, and 
is also more reliable than the mass of the p, since 
it can decay. In fig. [l]we show r /a as function of 
l//c sea for the three volumes. Even though part of 
this strong dependence of a on K sca could be reab- 
sorbed by improving the coupling || , it indicates 
the necessity to work at a fixed scale, whenever 
possible || . Moreover, while comparing 8 3 ■ 24 
and 12 3 • 24 there is a sizable difference at the 
lightest quark mass, the difference is completely 
negligible between 12 3 • 24 and 16 3 • 24. We can 
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Figure 1. ro/a as a function of 1/K sca . The solid 
errors are statistical and the dashed errors are 
systematic, as explained in 

summarise these results in a bound of L/rg > 3.2 
above which finite-size effects in the static quark 
potential are largely absent. Given the approxi- 
mate linear behaviour of Tq /a at the three lightest 
quark masses, we have attempted a chiral extrap- 
olation for the two largest volumes. In the chiral 
limit, we obtain a = 0.121(2) fm and a = 0.122(2) 
fm, for 12 3 • 24 and 16 3 • 24. 

One can rescale the linear Coulomb ansatz of 
the potential 



V(r) =V +ar 

r 



(2) 



e = tt/12, 
to eliminate a as follows 

[V(r)-V ]r = (1.65-e)(r/r -l)-e(r /r-l).(3) 

In fig. ||, we show the rescaled potential for 12 3 -24 
and for different values of K sca . The solid line 
doesn't represent a fit but simply a plot of eq. (||). 
In Full QCD, due to string breaking one expects 
the data to flatten out at large distances, below 
a value equal to twice the static meson mass. 
These limiting values have been calculated for 
K soa = 0.1390 (dashes) and K sca = 0.1395 (dots). 
With our present data, there is not conclusive ev- 
idence of string breaking up to r ~ 2.5r . A more 
chiral sea quark mass or a better way of extracting 
the potential are needed j|. At small distances, 
the data seems well described by the Coulombian 
part even if the data for the lightest quark masses 
tend to lie somewhat below the curve. This seems 
to favour a value of e > tt/12, as also observed in 
|?J for unimproved Wilson fermions. This is con- 
sistent with the influence of the dynamical quarks 
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Figure 2. Rescaled potential at V - 
the three lightest sea quark masses. 



12 3 ■ 24 for 



on the potential through the running of the strong 
coupling. We interpret this as a sign of unquench- 
ing and plan to carry out a quantitative analysis 
of the effect in a future publication, using the 
"correct" value of c sw |J . 

2. Light Spectrum 
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Table 2 

mps/my ratios for the different volumes. 

We have carried out a first complete analysis of 
the light hadron spectrum on the entire dataset, 
to gather complementary information to that ob- 
tained from the static quark potential. Ampli- 
tudes and masses have been obtained in a star- 
dard way by correlated least-x 2 fits to correlation 
functions. We have used throughout a double ex- 
ponential fitting function, taking into account the 
backward propagating state for mesons. We have 
found that the influence of the first excited state 
is not negligible as we approach the most chiral 
point. In order to estimate the excitation, we 
have fitted simulaneously the correlator "fuzzed" 
both at sink and source, denoted FF, and the cor- 
relator with local sink and source, denoted LL. 
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Figure 3. Meson masses for different volumes, in 
units of ro . 



The FF correlator allows a much faster isolation 
of the fundamental state than the LL correlator. 

In fig. H we show the vector and pseudoscalar 
mesons, as a function of 1/K sca , calculated at 
K vai _ K sea f Qr ^ e different volumes. To allow 
a significant comparison, we report all masses 
in units of ro, which reabsorbs the dependence 
of the lattice spacing on K sca . The plots show 
pronounced finite-size effects between 8 3 • 24 and 
12 3 • 24 as we move towards the most chiral point. 
On the other hand, between 12 3 • 24 and 16 3 • 24 
we find no significant discrepancy within statisti- 
cal accuracy at all values of the quark mass. A 
qualitatively similar picture is bourne out of the 
baryon sector. This behaviour confirms the one 
found for the static quark potential and we can 
extend to the mass spectrum the bound for finite- 
size effects of Lro <^ 3.2 already found there. 

One way to investigate sea quark effects in the 
spectrum is through the parameter J ||, which 
only requires the input of a physical V/PS mass 
ratio. In Full QCD, we need to fix K soa and vary 
the valence quark mass, interpreted as having 
strange flavour in a sea of light quarks. How- 
ever, our most chiral sea quark mass is only in the 
strange region. Moreover, the errors are amplified 
by the fitting process to extract the dmy/dmp S 
slope, especially at the lightest mass. In view 
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Figure 4. Vector vs pseudoscalar squared, in units 
of ro, for V = 12 3 • 24 at three values of K soa and 
comparison with quenched (K SGa = 0). 

of all this, it is not surprising that the values of 
J show no appreciable trend towards the exper- 
imental value as the sea quark mass decreases, 
as shown in r ef. [§. Another way to look for 
dynamical effects, equivalent in spirit but which 
bypasses any fitting procedure, is to look at the 
(mv,«ip S ) plane directly for a shift in the data, 
calculated at fixed K sca and different K val , as we 
approach the chiral limit. In fig. |] we show such a 
plot and also report the quenched result for com- 
parison (ft soa = 0), calculated at f3 = 5.7 and 
V = I2 3 • 24, which has a comparable lattice 
spacing. As the sea quark mass decreases, we do 
observe a significant, albeit small, trend towards 
the point (m^ , m^), i.e. the mesons whose flavour 
content resembles most closely that used in our 
simulation. Moreover, at the lightest sea quark 
mass, we can assert a significant shift compared 
to the quenched result. To quantify this quench- 
ing effect, we do need to keep the scale constant, 
an issue which will be addressed with the fully 
0(a) improved action ||. 
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